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ABSTRACT 

We have analyzed XMM-Newton Optical Monitor (OM) UV (180-400 nm) 
data for a sample of 33 galaxies. 30 are cluster member galaxies, and nine of 
these are central cluster galaxies (CCGs) in cooling flow clusters having mass 
deposition rates which span a range of 8 — 525 M© yr -1 . By comparing the ratio 
of UV to 2MASS J band fluxes, we find a significant UV excess in many, but 
not all, cooling flow CCGs, a finding consistent with the outcome of previous 
studies based on optical imaging data (McNamara &: O’Connell 1989; Cardiel, 
Gorgas, & Aragon-Salamanca 1998; Crawford et al. 1999). This UV excess is a 
direct indication of the presence of young massive stars, and therefore recent star 
formation, in these galaxies. Using the Starburst99 spectral energy distribution 
(SED) model of continuous star formation over a 900 Myr period, we derive star 
formation rates of 0.2 — 219 M© yr -1 for the cooling flow sample. For 2/3 of this 
sample it is possible to equate Chandra/XMM cooling flow mass deposition rates 
with UV inferred star formation rates, for a combination of starburst lifetime and 
IMF slope. This is a pilot study of the well populated XMM UV cluster archive 
and a more extensive follow up study is currently underway. 

Subject headings: cooling flows — galaxies:clusters:general — galaxies:elliptical and 
lenticular, cD — galaxies:stellar content — stars:formation — ultravioletrstars 
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1. Introduction 

The ’mystery’ of cooling flows (Fabian et al. 1984; Sarazin 1986; Fabian 1994) has baffled 
astronomers for over 20 years. Though the X-ray temperature often drops by a factor of three 
within the central 100 kpc (Allen et al. 2001), there is seldom enough mass at temperatures 
cooler than ~ 2 keV to be consistent with constant pressure cooling flow models. Recent 
XMM and Chandra results show that the distribution of temperature vs emission measure 
also does not match the prediction of these models (Peterson et al. 2003; Voigt & Fabian 
2004). While the XMM and Chandra results have reduced the implied cooling rate in most 
clusters there is still an apparent discrepancy of an order or magnitude or more between the 
X-ray cooling rates and the observed star formation rates (SFR) in many clusters. Contrary 
to conventional wisdom that there is no star formation at all in the central cluster galaxies, 
a number of contemporary studies have found significant star formation (Smith et al. 1997; 
Mittaz et al. 2001; McNamara et al. 2004a; O’Dea et al. 2004; McNamara 2004b). However 
we are still far from a complete reconciliation between X-ray mass deposition rates and 
alternatively derived star formation rates in the majority of cooling flow clusters. 

While many cluster cooling flows have been well observed in the radio and with emission 
line studies (Heckman et al. 1989, and others) there have been very few published studies 
of the UV properties of cooling flow clusters. In this study we use the readily available 
archival XMM Optical Monitor UV observations to derive the UV luminosities of a moderate 
sample of clusters. Since the dominant contribution to UV flux comes from short-lived main 
sequence stars, and the total number of these stars is directly proportional to the star 
formation rate (Kennicutt 1998) the UV one of the prime routes to understanding star 
formation. We attempt to address and quantify the connection between cooling flows and 
star formation in our pilot sample. Unless otherwise noted, this paper assumes a cosmology 
of H 0 = 70 km s -1 Mpc -1 , Um = 0.3, and Qa = 0.7 


2. Observations 

2.1. The XMM-Newton Optical Monitor 

The XMM Optical/UV Monitor Telescope is a modified 30cm Ritchey- Chretien tele- 
scope mounted alongside the X-ray mirror modules (Mason et al. 2001) which allows simul- 
taneous X-ray and UV/optical observations of XMM targets within the central 17 square 
arcminute region of the X-ray field of view. The 256 x 256 pixels, have a ~ 1" . The sen- 
sitivity of the Optical Monitor (OM) is comparable to a ground based 4m telescope. This 
study primarily utilizes the UVW1 filter (220-400 nm), and the UVW2 filter (180-260 nm), 
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and the UVM2 filter (200-280 nm) for calibration purposes. Previous XMM OM studies of 
clusters have been limited to the detection of UV mission in Abell 1795 (Mittaz et al. 2001). 


2.2. Galaxy Sample 

Our pilot sample consists of 33 galaxies and eleven clusters. Three of them (M49, M60, 
and M86) were selected because they possess IUE photometric information (Burstein et al. 
1988) and were used to compare OM to IUE fluxes. 21 galaxies in the IC 1860 and Coma 
fields were used to calibrate the relationship between UV and IR luminosities in passively 
evolving early type and cD galaxies. These two clusters were chosen because they had 
archival UVW1 data, are at low redshift, have numerous early type galaxy detections with 
2MASS survey data, and do not have cooling flows. We utilized two UVW1 observations of 
Coma in our study, one centered on the cD galaxies, and one which is off-center and contains 
numerous early-type galaxy detections within the viewing field. 

The remaining nine clusters in our sample were chosen because of significant X-ray de- 
termined mass deposition rates, general overlap with the samples of McNamara & O’Connell 
(1989) and Cardiel, Gorgas, & Aragon-Salamanca (1998), moderate to low redshift, inclu- 
sion in the 2MASS galaxy sample, and having archival XMM UVW1 data. The majority of 
clusters in our sample are low redshift objects (z < 0.1), but we have also included three 
moderate redshift clusters (z ~ 0.25) which are known to possess very strong cooling flows. 

Figure 1 shows one of the cooling flow clusters in our sample, with XMM-Newton UVWl 
contours overlayed. Table 1 lists the objects in our current sample, their redshifts, relevant 
filters, and OM exposure times. 


3. Photometry 
3.1. Data Analysis 

We retrieved pipeline-processed UV data from the XMM-Newton Science Archive and 
combined observations with the same pointing. Flux determinations for our calibration 
sample used a circular aperture of 7" radius, a good approximation of the IUE aperture 
(Burstein et al. 1988), and the size of the fixed 2MASS photometric aperture. Photon 
counts in the aperture were corrected for background using a 35" radius background region 
taken from an outlying area of the chip in the same observation. Final fluxes were determined 
by employing the OM counts-to-flux conversions, and errors were assessed via \/N photon 



-4- 


* ^ 


statistics. 

The results from UVW2 and UVM2 observations of M49, M60, and M86 were then 
compared with the IUE determined flux value of each object (Burstein et al. 1988) for 
the bandpass of the appropriate filter. Our flux measurements are consistent with theirs, 
indicating that the Optical Monitor does not have a red leak, and that our method of 
determining UV fluxes is reliable. 

We then determined the UVW1 fluxes for the remainder of our sample (see Table 1), 
for a 7" radius region in the method described above. The UVW1 data were corrected for 
average galactic extinctions, Cardelli, Clayton, & Mathis (1989) and J band extinction values 
were retrieved from the NASA/IPAC Extragalactic Database (NED). 


3.2. Calibration 

The low redshift calibration sample has 23 galaxies in IC 1860 and Coma with 2MASS 
J band photometric measurements. Two of these galaxies exhibited high UV fluxes with 
respect to the rest of the calibration sample, and were described as star forming in a FOCA 
study Donas, Millard, & Laget (1995) and thus are not appropriate for inclusion in our 
UV/IR luminosity calibration of passively evolving early type galaxies, and were removed 
from this study. 

A least squares fit between the UVW1 and J band luminosities for the remaining 21 
early type galaxies be is well described by 


Luv = Ci x Lj^ 2 (1) 

Where Luv is the UVW1 luminosity and Lir is the 2MASS J band luminosity, inside 
7" radius apertures with C\ = 0.090 ± 0.007 and C 2 = 0.708 ± 0.007 respectively (Figure 
2). We use this correlation to predict the UV light from non-starforming, passive galaxies 
assuming that the UV light from such objects is produced by old evolved stars coming from 
the same population that produces the IR light. 


3.3. Photometric Results 

The UV luminosity excesses of the cD galaxies in our nine cooling flow clusters were 
calculated by subtracting the expected luminosity (based on the fit obtained in Section 3.2) 
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from the measured value and Figure 3 shows the J band luminosity of each galaxy in our 
sample plotted against its UV/IR luminosity ratio. We have not attempted to estimate and 
correct for internal dust absorption, and thus these excesses provide a lower limit on the UV 
emission in these clusters. 

Starburst99 (Leitherer 1999) redshifted models corresponding to continuous star forma- 
tion at 1 M© per year over a 900 Myr period, using an IMF powerlaw index of 3.3, were 
folded through the UVW1 filter to produce the UVW1 luminosity. The 900 Myr model was 
chosen is closest to the average lifetime of a cooling flow (Allen et al. 2001). Based on the 
relationship between UV luminosity and star formation rate (Kennicutt 1998), our measured 
UV luminosity excesses for each cluster were then used to estimate the per-year star forma- 
tion rate of each galaxy over 900 Myr. 2, and compared to X-ray derived mass deposition 
rates from Voigt & Fabian (2004), Peterson et al. (2003), and R. Mushotzky (2005, private 
communication). Also in Table 2 are estimates of the duration of continuous star formation, 
and best fitting IMF powerlaw index. These values were derived by setting the UV mass 
deposition rate equal to the X-ray mass deposition rate and determining which model best 
described that scenario. 


4. Surface Brightness 

We fitted the flux in concentric circles around the emission peak of each galaxy and to 
a radial surface brightness profiles of the form: 

/ „2\ 1/2— 3/3 

I(r)= io ( 1 + 3J + Ib (2) 

Where 7 0 is the normalization, ro is the core radius, and Ib is a constant representing the 
intensity contribution of the background. 

The surface brightness profiles are well fit by these models, cf. Figure 4. The fit 
parameters are in Table 3. 


5. Summary and Discussion 

We find a significant UV luminosity excess in the central galaxies of most cooling flow 
clusters. Which implies that star formation has recently been or is currently occurring 
most likely as a direct result of radiative cooling of the intracluster medium (ICM). These 
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conclusions are qualitatively consistent with the findings of previous optical imaging studies 
of star formation in cooling flow clusters (McNamara & O’Connell 1989; Cardiel, Gorgas, & 
Aragon-Salamanca 1998; Crawford et al. 1999). 

In the absence of multiwavelength spectral information, it is difficult to constrain a 
particular model of the star formation history of these galaxies. However we note that for 
2/3 of our sample it is possible to find some combination of the parameters of IMF slope and 
starburst lifetime which allow us to equate X-ray determined mass deposition rates with UV 
inferred star formation rates. 

In none of the OM, Chandra or HST images is there evidence of a significant excess 
central flux, therefore we expect the AGN UV contribution to our signal is small. Since the 
calibration sample consists entirely of low redshift early type galaxies, there is a possible 
systematic error into our calculations of the SFRs of the three moderate redshift (z ~ 0.25) 
clusters in our sample (Brown et al. 2003). In future work we will calibrate a higher redshift 
passive E galaxy sample. 

While there is significant UV flux in 2/3 of our cooling flow sample, 1/3 of the objects 
with significant cooling rates from modern data do not have detectable UV emission, and 
the ratio of UV light to X-ray cooling rate (Figure 5) shows a wide range. The origins of 
this scatter are unclear and a larger UV sample, with direct comparison of X-ray and radio 
imaging and IR data, is necessary to get a clearer picture of the origin of the UV light. The 
discovery of copious UV emission from several cooling flow clusters confirms the previous 
work and shows that a significant fraction of the cooling gas in clusters can indeed form 
stars. 

This paper is a pilot study and we are currently undertaking an expansion of this 
project toward higher redshift, utilizing more of the OM UV filters and including the many 
additional clusters that remain in the archive. Through this expansion we hope to achieve a 
more thorough characterization of recent star formation in cooling flow clusters of galaxies. 


This work was made possible by a NASA Graduate Student Research Program (GSRP) 
fellowship, NGT5-140. 

This research has made use of the NASA/IPAC Extragalactic Database (NED) which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Administration. 
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Fig. 1. — HST WFPC2 1780 second exposure image of Abell 1795. This image has been 
overlaid with UV contours from an XMM-Newton Optical Monitor UVWl filter (220-400 
nm) observation. Logarithmic contours in black were produced with a 2499 second exposure 
image, and have count values of 14. 25, 44.75, and 80. The large splotch of UV contours to 
the lower right of the cD galaxy is due to scattered light. 
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Fig. 2. — UV vs. IR luminosity for non-cooling flow cluster galaxies in our sample. The line 
indicates the best fitting relationship between UV and IR luminosity for passively evolving 
and non-cooling flow cD galaxies, represented by squares and triangles, respectively. 
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J Band Luminosity vs. UV/IR Ratio 



Fig. 3. — UV/IR luminosity ratio of selected galaxies vs. their 2MASS J band luminosities. 
All luminosities were calculated in 7 arcsecond radius regions, and corrected for galactic 
extinction. The small squares denote field ellipticals (M49, M60, and M86), while the larger 
squares indicate passively evolving early-type cluster member galaxies at low redshift. The 
triangles represent low redshift non-cooling flow cD galaxies, and the diamonds indicate 
cooling flow cluster cDs, all of which are at low redshift (z < 0.1) except for the three points 
farthest to the right, which have redshifts of ~ 0.25. The high UV/IR ratio of the cooling 
flow clusters is a strong indication of the existence of hot young stars and thus a direct 
indication of copious recent star formation. 
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Fig. 4. — Surface Brightness Profile. Radial surface brightness profile of Abell 1795 
made from UVW1 (220-400 nm) data accumulated in ~ 1" annular bins. A solid line traces 
the best (3 model fit. The horizontal dotted line represents the best fitting background value. 
The best fitting (3 model has a core radius of 4.09 " and beta value of 0.703 
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Excess UV luminosity vs. til 
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Fig. 5. — Excess UV luminosity vs. M. Luminosity excesses were determined by subtracting 
the best fit ratio of UV light due to old stars (Figure 2) from the measured UV luminosity 
of each cluster. X-ray mass deposition rates were taken from recent papers (Peterson et al. 
2003; Voigt & Fabian 2004). 



Table 1. XMM-Newton Optical Monitor Observations 


- Object 

z 

Filter 

Exposure [s] 

IC 1860 

0.0229 

UVW1 

5000 

PKS 0745-191 

0.1028 

UVW1 

3999 

Hydra A 

0.0538 

UVW1 

3000 

Zw3146 

0.2906 

UVW1 

4000 

NGC 4406 (M86) 

-0.0008 

UVM2 

3500 



UVW1 

3600 

NGC 4472 (M49) 

0.0033 

UVM2 

1000 



UVW1 

2000 

NGC 4649 (M60) 

0.0037 

UVW2 

2900 



UVW1 

2499 

Coma (center) 

0.0231 

UVW1 

3900 

Coma (outskirts) 


UVW1 

1498 

Abell 1795 

0.0633 

UVW1 

2499 

Abell 1835 

0.2520 

UVW1 

2000 

MS 1455.0+2232 

0.2578 

UVW1 

2000 

Abell 2029 

0.0779 

UVW1 

1200 

Abell 2052 

0.0345 

UVW1 

2000 

MKW3s 

0.0453 

UVW1 

2002 
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Table 2. 


Cluster 

z 

UV Excess 
[10 43 erg/s] 

M (900 Myr) * 
[M e /yr] 

M (X-ray) 
[M 0 /yr] 

T s i 

[Myr] 

a 

PKS 0745-191 

0.1028 

11.94 ±a7 

219 ±12 

317±|| b 

200 

3.3 

Hydra A 

0.05384 

0.88 ±002 

16.9 ±0 - 4 

8 « b 

900+ 

2.35 

Zw3146 

0.2906 

8.47 ±0 - 4 

121 ±6 

525lgo 

11 

3.3 

Abell 1795 

0.06326 

O.85 ±002 

16.1 ±0 - 4 

30 + l h 

100 

3.3 

Abell 1835 

0.251982 

11.41 ±0 - 4 

174 ±7 

34lg b 

900 + 

2.35 

MS 1455.0+2232 

0.2578 

3.36 ±a4 

51 ±6 

200i?g 

150 

3.3 

Abell 2029 

0.07795 

0 i8±° °5 

3 4±°- 9 

30 t\l b 

7 

3.3 

Abell 2052 

0.03446 

0 i9±°.°°7 

0 4±°.i 

iool 2 ° c 


. . . 

MKW3s 

0.04531 

O 

O 

H- 

P 

o 

t— * 

0 2±°- 2 

c 




a Using an IMF powerlaw of 3.3. Star formation rates for a Saltpeter IMF fall close to 
half of this value. 

b Voigt & Fabian (2004). 

c Peterson et al. (2003). 
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Table 3. Beta-Model Fits 


Object 

Rcore ["] 

P 

V 

/ B a 

Reduced x 2 

PKS 0745-191 

2 83± 0 04 

0.766 ±oon 

625. 3 ±124 

481.0 ±1 - 6 

0.691 

Hydra A 

3.449 ±0 °oi 

O.8387 ±00002 

4829. 1 ±24 

356. 59 ±006 

4.200 

Zw3146 

6.2l ±001 

2 o54 ± 0 007 

1119. 2 ±4 3 

413. 2 ±04 

0.947 

Abell 1795 

4 09±° °i 

0 703± 0 °°2 

2334. 8 ±13 ' 4 

336. 6 ±0 ' 3 

0.946 

Abell 1835 

3.65 ±a41 

1.253 ±0165 

2880.7 ±1530 

349. 6 ±L8 

0.514 

MS1455. 0+2232 

4 99±°.°3 

1.006 ±oo ° 7 

446. 9 ±35 

337.8 ±0 ' 4 

0.657 

Abell 2029 

3.68 ±002 

O.499 ±ao01 

613. 2 ±3 ' 5 

449. 7 ±a6 

1.101 

Abell 2052 

1 68±°-28 

0 402±° 022 

1593. 1 ±I49 ' 6 

426. 2 ±74 

1.912 

MKW3s 

1.28 ±0 ' 33 

O.352 ±0022 

1018. 0 ±138 9 

398. 6 ±71 

1.110 


a Surface brightness / in units of 10 17 ergs sec 1 cm 2 arcsec 



